Isospin Chemical Potential and the QCD Phase Diagram at Nonzero Temperature 

and Baryon Chemical Potential 
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We use the Nambu-Jona-Lasinio model to study the effects of the isospin chemical potential 
on the QCD phase diagram at nonzero temperature and baryon chemical potential. We find that 
the phase diagram is qualitatively altered by a small isospin chemical potential. There are two first 
order phase transitions that end in two critical endpoints, and there are two crossovers at low baryon 
chemical potential. These results have important consequences for systems where both baryon and 
isospin chemical potentials are nonzero, such as heavy ion collision experiments. Our results are in 
complete agreement with those recently obtained in a Random Matrix Model. 



Introduction. - QCD at nonzero temperature and 
baryon chemical potential plays a fundamental role in 
many different physical systems. Two important ones are 
neutron stars, which probe the low temperature and high 
baryon chemical potential domain, and heavy ion colli- 
sion experiments, which probe the high temperature and 
low baryon chemical potential domain. In the last few 
years, various models have been used to predict the main 
characteristics of the QCD phase diagram at nonzero 
temperature and baryon chemical potential Q . The exis- 
tence of color superconducting phases at low temperature 
and high baryon chemical potential 13 j as well as the 
presence of a tricritical point at intermediate tempera- 
ture and baryon chemical potential [5| are among the 
most important results. These features have numerous 
phenomenological consequences. Some of the effects of a 
nonzero isospin chemical potential have also been stud- 
ied, but only in the low temperature and high baryon 
chemical potential domain [1, 0, H, OjJ ■ However 
systems such as heavy ion collision experiments that ex- 
plore the high temperature and low baryon chemical po- 
tential domain also have a nonzero isospin chemical po- 
tential. Therefore, there is a clear need to study the ef- 
fects of a nonzero isospin chemical potential on the whole 
QCD phase diagram at nonzero temperature and baryon 
chemical potential. 

Most of our knowledge of the QCD phase diagram at 
nonzero temperature T is restricted to either zero baryon 
chemical potential /i^, or to zero isospin chemical po- 
tential /i/. At /iB=/i/=0 numerical lattice simulations 
and effective theories predict that the ground state cor- 
responds to a hadronic phase at low temperature and to 
a quark-gluon-plasma phase at high temperature. For 
nonzero quark masses, there is no order parameter that 
distinguishes between these two phases. For QCD with 
two flavors, a crossover is expected at T ~ 170 MeV [T^ . 
This crossover extends into the phase diagram at nonzero 
baryon and isospin chemical potentials. 

Numerical lattice simulations have explored the high 
temperature and small baryon chemical potential domain 
lis, 14, 15, 16, 17, 18]. It was found that these high 
temperature crossover lines at /ib Aqcd and = 0, 
and at /is = and fij much smaller than the pion mass 



are identical 0, . 

At zero baryon chemical potential, both effective the- 
ories and numerical lattice simulations predict the exis- 
tence of a superfluid pion condensation phase for high 
enough a*/ El IH 113, M IH IH H . At zero tempera- 
ture a second order phase transition at a critical isospin 
chemical potential, /x"'*, equal to half the pion mass sep- 
arates the hadronic phase from the pion condensation 
phase. When the temperature is increased, this second 
order phase transition line ends in a tricritic al p oint and 
the phase transition becomes first order |23.l23 |. 

At nonzero baryon chemical potential, standard nu- 
merical lattice simulations do not work. Therefore our 
knowledge of the QCD phase diagram at nonzero T 
and /is relies exclusively on effective theories, such as 
Nambu-Jona-Lasinio and Random Matrix models. At 
temperatures smaller than a few tens of MeV, an in- 
crease in fj,B leads to a crystalline LOFF phase 0- If 
/is is further increased, the ground state corresponds to 
a color superconductor. Both of these phase transitions 
are of first order. If the temperature is increased to a few 
tens of MeV the LOFF and color superconducting phases 
disappear, and a first order phase transition directly sep- 
arates the hadronic phase from the quark-gluon-plasma 
phase. At zero quark mass, this first order line ends in 
a tricritical endpoint at T ~ lOOMeV, where a second 
order phase transition starts 0, |^ . 

Finally, the only available theoretical study of the 
whole phase diagram at nonzero T, fj,B, and u j has been 
performed using a Random Matrix model |24| . It was 
found that a small isospin chemical potential induces two 
first order phase transitions at low T that end in two crit- 
ical endpoints, and there are two crossovers at low /i^. 
Because of the phenomenological implications of these re- 
sults, it is essential to try to reproduce them within other 
models. 

In this letter, we study the QCD phase diagram at 
nonzero T, /is, and /i/ for two quark fiavors of equal 
mass ni in the Nambu-Jona-Lasinio model. As a first 
step and for simplicity, we shall not study the domain 
of low temperature and high baryon chemical potential 
where the ground state corresponds to a LOFF crystal or 
to a color superconductor. We shall also restrict ourselves 
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to small isospin chemical potential, ^/ < /zj''''. A more 
complete analysis will be published elsewhere. 

Phase Diagram in the Nambu-Jona-Lasinio Model. - 
A crucial observation made in is that when both 
fiB 7^ and ^/ 7^ in QCD, there is no reason to expect 
that the quark-antiquark condensates are equal for each 
flavor. Indeed if we define fiB = ^(Mu + Md) and /i/ = 
— fid), the QCD Lagrangian can be written as 



^QCD = ^ tl^f {ip - m + fifjo)ipf, 



(1) 



f=u,d 



and there is no symmetry in the QCD Lagrangian that 
constrains {uu) to be equal to {dd). Therefore, we have to 
consider the quark-antiquark condensates for each flavor 
separately. In this letter, we shall concentrate on three 
observables: The quark-antiquark condensates for each 
flavor cr„ = {uu) and <Td = {dd), as well as the pion 
condensate p = ^{u^^d — d'y^u). 

We use the same Nambu-Jona-Lasinio model with 
an instanton-induced four-fermion interaction as Berges 
and Rajagopal who studied the QCD phase diagram at 
nonzero T and fis in 5] • After the standard introduction 
of bosonic flelds via a Hubbard-Stratonovich transforma- 
tion and integration over the fermion fields |^, |25| , we 
find that the mean-field free energy is given by 



-Tr log 



< + 2p' 
hu 



(2) 



where 



hf = {iujri + M/)7o + iiry + m + F^{p)af, (3) 

with ujn = {2n + 1)ttT, and the form factor 

F{p)=Ay{f + A') (4) 

is introduced to mimic the effects of asymptotic freedom 
In the free energy 10), we have only kept the poten- 
tially non- vanishing condensates cr„, ad, and p. We follow 
Berges and Rajagopal and take the scale A = 0.8 GeV 
and the coupling constant Gi = 6.47/A^ which are rea- 
sonable phenomenological choices j^. They correspond 
to Ou = cfd = 0.4 GeV at m=T=fiB=P'i=^- 

In order to proceed we have to compute the excita- 
tion energies in ((SJ, and in general solve a fourth order 
equation in uJn- We have not been able to find a short 
expression for the general solution. We shall therefore an- 
alyze the most relevant particular cases separately. First 
at Mb = (i.e. fiu = —fJ-d), relying on lattice simulations 
[l8j , we can assume that cr„ = (Xd = a. We find that the 
excitation energies are given by 



E± = J{E±fiir + F^ip)p^ 



(5) 



where E = y/jp + (m + F'^{p)(t)^. Thus the free energy 
(O becomes 
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^ -10 



dp p 



E±+2Tlog(l + e-^^/^ 



(6) 



This is exactly the same free energy as the one that has 
been studied for the phase diagram of QCD with two 
colors in l^^. It leads to a Bose- Einstein condensation 
phase where p ^ when pj is larger than half the pion 
mass |2^. Therefore at pB = 0, the Nambu-Jona-Lasinio 
model agrees with the results obtained in lattice simula- 
tions and in effective theories 0, I^O, 0, 0] • 

We then study pB ^ and pi < pf^^, which corre- 
sponds to the most relevant phenomenological situations, 
i.e. when the pion condensate vanishes. The free energy 
(O then separates into a sum over each flavor 



E 



(7) 



■-^ / dpp^ 



Ei 



2Tlog(l + e-^±/^) 



where 



E: 



f _ 



+ (m + F'^{p)afY ± pf. 



(8) 



The free energy {TJ has two remarkable properties. It is 
even in pu and pd separately and it can be expressed as a 
sum over the different quark flavors. Both of these prop- 
erties are also found in the Random Matrix model stud- 
ied in |23| and are essentially responsible for the striking 
changes in the phase diagram. 

The evenness of the free energy in pu and pd implies 
that the crossover line that separates the hadronic phase 
from the quark-gluon-plasma phase in the pB — T plane 
at /i/ = coincides with the corresponding crossover line 
in the pi — T plane at /is = 0. This property was indeed 
found in numerical lattice simulations 15, 18]. 

Since Vl — J2f=u d^fiPf)^ the free energy is mini- 
mized by minimizing each f2 j separately. But each j is 
equal, up to a factor two, to the free energy studied _by 
Berges and Rajagopal for zero diquark condensate in |5(. 
Therefore, for each flavor separately, the domains where 
aj ^ 0.4 GeV and aj << 0.4 GeV are separated by a 
flrst order line that ends in a critical endpoint where a 
crossover lines starts. Now since flu{pu) = ^u{pb + Pi) 
and fld{pd) — ^u{pb — Pi) , the whole phase diagram for 
QCD with two flavors in the pB — T plane at nonzero pi 
corresponds to a superposition of two of the usual phase 
diagrams at pi = shifted by 2pi. This is illustrated 
in Figure 1 for to = 10 MeV and for both pj = and 
Pj = 30 MeV. There are now two first order lines that 
start at r = and which end in a critical endpoint at 
T ^ 65 MeV. The temperature of the critical endpoint 
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FIG. 1: Phase diagram in the /is-T-plane for a quark mass m = 10 MeV and an isospin chemical potential /ij = 0, and 
fii = 30 MeV, respectively. At low temperature a first order phase transition takes place at the full line that ends in the critical 
endpoint. The dotted curves depict the crossover behavior. The condensates that are not displayed are of the order of the 
quark mass. This phase diagram can be trusted only above temperatures of a few tens of MeV since we did not consider the 
possibility of a crystalline LOFF phase or a color superconductor. The temperature of the critical endpoint is not affected by 
the isospin chemical potential. 



is not affected by the isospin chemical potential. Two 
crossovers emerge from the critical endpoints and inter- 
sect on the fiB — axis at T ^ 200 MeV. Notice that 
the very low temperature part of this phase diagram can- 
not be trusted since we did not consider the possibility 
of a crystalline LOFF phase or a color superconductor. 
Nevertheless the phase diagram above a few tens of MeV 
is correct since the LOFF and the color superconducting 
phases disappear for such temperatures. 

Conclusions and Discussion. ~ In this letter we have 
used a Nambu-Jona-Lasinio model to show that the in- 
troduction of an isospin chemical potential leads to qual- 
itative changes in the QCD phase diagram at nonzero 
temperature and baryon chemical potential. 

First, in agreement with lattice simulations [Tsl IT^ 
and the Random Matrix model analyzed in we find 
that the crossover line at low hb and /i/ = is identical 
to the crossover line at low iij and iib = 0- 

Second, at low temperature, there are two first order 
phase transitions that end in two critical endpoints, and 
there are two crossovers at low baryon chemical poten- 
tial. All these results are in complete agreement with the 
Random Matrix model studied in ^2?\. The new lattice 
techniques developed to study the small /is behavior at 
m — can be used to study this pair of crossover lines. 

The existence of this pair of crossover lines is impor- 
tant for heavy ion collision physics. First, as was shown 
in [23, m, I23I, the critical endpoints have definite sig- 
natures that can be observed in these experiments. The 



effect of the isospin chemical potential is twofold: It dou- 
bles the number of critical endpoints and pushes one of 
them to lower /i^. The second effect makes the pres- 
ence of the critical endpoint easier to see in heavy ion 
collision experiments. Second, the transition from the 
quark-gluon-plasma to the hadronic phase will be softer 
at nonzero fij than at fii = 0, since the system has to go 
through two crossover lines instead of one. 

Finally if the strange quark is included, the hadronic 
phase and the quark-gluon-plasma phase are expected 
to be separated either by a crossover or by a first order 
phase transition at zero fj,i, depending on the precise 
value of the strange quark mass ^3 ■ The effect of a 
small strange quark mass is to push the critical endpoint 
towards higher T. We expect that the introduction of 
a nonzero fij at a physical value of the strange quark 
mass will generate two of these separation lines: Either 
two crossovers or two first order phase transitions will 
be present at high temperature and low baryon chemical 
potential. 
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